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The  authors  present  their  experimental  and  theoretical  results  on  nonlinear  absorption  of  light  in 
InAs.  The  nonlinear  variation  of  output  intensity  as  a  function  of  input  intensity  and  time  are 
calculated  by  solving  four  coupled  rate  equations  simultaneously.  Ah  required  quantities,  including 
two-photon  absorption,  free-carrier  absorption.  Auger  and  radiative  recombination  lifetimes,  and 
intrinsic  carrier  densities,  have  been  obtained  from  the  underlying  band  structures.  The  calculated 
thickness  and  energy-dependent  output  intensities  in  InAs  agree  very  well  with  the  values  measured 
in  their  pump-prohe  experiments.  ©  2006  American  Institute  of  Physics.  [DOT  10.1063/1.2363970] 


Nonlinear  absorption  in  semiconductors  continues  to  be 
of  scientihc  and  technological  interest.^  Technologically, 
nonlinear  absorption  (NLA)  is  being  used  in  a  number  of 
applications,  including  frequency  conversion,  optical 
switches,  and  optical  limiting.*  The  importance  of  free- 
carrier  effects  in  controlling  the  infrared  nonlinear  optical 
properties  has  long  been  realized. Consequently,  consider¬ 
able  efforts  were  undertaken  to  understand  the  NLA  and 
evaluate  the  parameters  such  as  two-photon  absorption 
(TPA)  coefficient,  free-carrier  absorption  (FCA)  coefficient, 
refraction  cross  section,  and  lifetimes  in  semiconductors.^  *** 
However,  extracting  parameters  from  the  measured  values 
usually  requires  careful  design  of  experiments  to  isolate  and 
remove  the  effect  of  other  parameters  on  the  observed  re¬ 
sults.  Even  when  evaluated  at  appropriate  limits,  the  extrac¬ 
tion  ignores  the  complicated  interdependence  of  these  pa¬ 
rameters  on  temperature  and  carrier  density.  For  example, 
TPA  and  FCA  depend  on  carrier  density,  which  in  turn  de¬ 
pends  on  temperature.  Consequently,  the  extracted  values 
differ  substantially  from  each  other.'*^®’^’***  In  addition,  it  is 
not  clear  whether  the  values  extracted  in  certain  limits  can  be 
used  in  modeling  of  high  intensity  light  propagation  where 
the  photoexcited  carrier  distribution  is  far  from  equilibrium. 
There  have  also  been  several  intuitive  calculations  to  evalu¬ 
ate  these  parameters  using  either  effective  mass  bands  and 
k-p  bands  with  adjustable  matrix  elements.  ’  The  trends 
are  often  obtained  correctly.  However,  for  complete  under¬ 
standing  and  modeling  of  light  propagation  with  improved 
predictability,  an  accurate  evaluation  of  these  parameters 
from  full  band  structures  needed  to  be  carried  out. 

In  this  letter,  we  report  the  results  of  transmission  mea¬ 
surements  carried  out  with  (a)  pump  only  and  (b)  the  pump- 
probe  technique  on  a  well-characterized  InAs  sample  along 
with  rigorous  calculation  of  time  (t)  and  space  (z)  dependent 
transmitted  pump  and  probe  intensities.  All  physical  quanti¬ 
ties  including  TPA  and  FCA  coefficients,  Auger  recombina¬ 
tion  (AR),  and  radiative  recombination  (RR)  lifetimes  are 
calculated  as  functions  of  temperature  (T)  and  carrier  density 
(N)  using  full  band  structures.  The  calculated  value  of  the 
transmitted  probe  intensity  at  the  exit  surface  explains  the 
experimentally  observed  trend  very  well  and  an  excellent 


quantitative  agreement  is  obtained  when  a  Shockley-Read- 
Hall  (SRH)  mechanism  is  also  included. 

In  our  pump-probe  experiment,  the  charge  carriers  were 
generated  using  TPA.  In  addition  to  the  convenient  availabil¬ 
ity  of  below  band  gap  lasers,  the  TPA  used  here  ensures  that 
the  entire  thickness  of  the  sample  is  utilized  for  carrier  gen¬ 
eration  at  low  incident  irradiance.  The  4.8  fim  pump  source 
was  obtained  by  frequency  doubling  a  transversely  excited 
atmospheric  CO2  laser  to  provide  128  ns  (full  width  at  e“*  of 
the  maximum)  duration  pulses  with  energies  up  to  10  mJ. 
First,  we  carried  out  a  set  of  NLA  experiments  at  room  tem¬ 
perature  with  only  the  pump  present. 

In  Fig.  1,  the  intensity  transmission  of  the  sample  is 
plotted  as  a  function  of  incident  intensity  for  four  different 
InAs  sample  thicknesses  (<7=0.11,  0.79,  1.21,  and  1.93  mm). 
With  increasing  intensity,  more  carriers  are  created  by  TPA, 
and  the  resulting  FCA  decreases  the  output  intensity  nonlin- 
early  with  input  intensity.  As  the  thickness  of  the  sample  is 
increased,  the  path  length  of  the  light  is  longer  and  conse¬ 
quently  more  NLA  takes  place,  as  shown  in  Fig.  1. 

The  decay  time  of  the  TPA  excited  carriers  is  studied 
using  a  pump-probe  method.  The  pump  is  the  4.8  fim  beam 
described  above  and  the  probe  is  the  cw  laser  beam  obtained 
by  a  grating  tuned  CO  laser  operating  at  5.3  fim  with  power 


FIG.  1 .  Compari.son  of  tran.smitte<l  pump  intensity  measuretl  (data  points)  at 
300  K  in  10*^  cm“^  n-doped  InAs  with  that  calculated  (solid  lines)  for  four 
thicknesses. 
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FIG.  2.  Comparison  of  transmission  intensity  measured  (data  points)  at 
300  K  in  10*^  cm“^  «-doped  InAs  with  that  calcutated  (solid  lines)  for  three 
beam  energies. 


of  100  mW.  The  probe  beam  was  chopped  by  a  mechanical 
shutter  to  a  5  ms  duration  pulse.  A  dichroic  filter  (coated  for 
high  reflection  at  4.8  /ttm  and  antireflection  at  5.3  fim  for 
45°  angle  of  incidence)  was  used  to  combine  the  pump  and 
probe  beams  which  were  focused  on  the  sample  with  the 
same  200  mm  focal  length  lens  to  spot  sizes  of  378  and 
166  /u-m,  respectively  (beam  radius  at  e"'  of  maximum).  The 
sample  studied  is  a  1-mm-thick  InAs  wafer  with  antireflec¬ 
tion  coating.  The  probe  transmission  was  measured  using  a 
HgZnCdTe  detector  with  a  rise  time  less  than  1  ns,  and  the 
incident  and  transmitted  pump  energies  were  measured  using 
pyroelectric  detectors. 

Figure  2  displays  the  time  dependence  of  the  transmitted 
probe  intensity  (data  points)  through  the  sample  for  three 
different  pump  beam  energies.  At  the  onset  of  the  pump 
pulse  (r=0),  TPA  in  the  material  generates  electron-hole 
pairs.  Thus  created  free  carriers  absorb  probe  beam,  by  FCA, 
and  cause  a  decrease  in  probe  transmission.  The  transmission 
is  minimum  at  the  peak  of  the  pump  pulse.  As  the  free  car¬ 
riers  recombine,  their  density  decreases  and  the  FCA  de¬ 
creases.  Consequently,  the  probe  transmission  increases  to 
eventually  reach  its  initial  value.  With  an  increase  in  the 
pump  energy,  from  25  to  124  fiJ,  more  carriers  are  created  at 
the  end  of  the  pulse  and,  consequently,  the  minimum  in 
probe  transmission  is  reduced  further. 

For  a  detailed  understanding  of  the  mechanisms  affect¬ 
ing  the  light  propagation  in  the  material,  a  complete  solution 
to  the  following  rate  equations  is  needed. 


cients,  respectively,  and  cr  (cr*)  denote  the  sum  of  (difference 
in)  the  cross  sections  for  FCA  arising  from  phonon  emission 
and  absorption.  The  total  recombination  rate  is  the  sum 
of  the  AR  rate  (t^r),  RR  rate  (trr),  and  SRH  rate  (t^rh)- 
is  the  intrinsic  carrier  density  at  a  given  T.  We  have  consid¬ 
ered  T  increase  due  to  the  FCA  and  TPA,  which  has  not  been 
included  in  previous  studies.  We  find  that  the  increase  in  T 
can  be  significant  (as  much  as  150  K  for  moderate  inten- 
sites).  It  is  important  to  realize  that  Nq,  a,  f3,  cr,  and  Tr  are  T 
dependent.  In  addition,  (3  and  Tr  are  dependent  on  N,  which 
in  turn  depends  on  T.  Moreover,  the  carrier  density  created 
by  the  pump  beam  affects  the  absorption  of  the  probe  beam. 
Hence,  Eqs.  (l)-(4)  have  to  be  solved  self-consistently.  For 
the  chosen  value  of  incident  pump  intensity  Ip  {t=Q,r,z 
=  0)  at  a  radial  distance  r,  the  time-dependent  probe  intensity 
/pj  it,r,z)  is  evaluated  at  the  exit  side  iz=L)  of  the  sample  to 
compare  with  experiments.  Since  the  photon  energy,  hv,  is 
smaller  than  the  band  gap  of  the  materials  considered,  a  is 
zero.  The  self-consistent  solution  requires  that  Nq,  (3,  cr,  and 
Tr  be  calculated  at  all  T  and  N  values  of  interest.  We  calcu- 


of  10'^  cm  10' 


and 


late  them  in  the  range 
80  K=sT«500  K. 

The  TPA  and  FCA  coefficients,  RR  and  AR  lifetimes, 
Fermi  functions,  and  intrinsic  carrier  concentrations  are  cal¬ 
culated  using  a  hybrid  pseudopotential  tight-binding  Hamil¬ 
tonian  that  includes  all  long-range  interactions."  In  general, 
(3  is  assumed  to  be  N  and  T  independent.  However,  we  see 
that  the  Moss-Burstein  effect  offers  serious  limitations  to  p 
at  moderate  concentrations  and  low  temperatures.  For  a 
given  value  of  (equal  electron  and  hole  density)  N  and  T,  we 
first  obtain  the  quasi-Fermi  level,  as  appropriate  for  the  non¬ 
equilibrium  studies  considered  here.  Then,  using  appropriate 
electron  and  hole  distributions  (for  final  and  initial  states, 
respectively),  (3  is  calculated."  We  find  that  (3  (at  4.8  /i,m) 
remains  nearly  constant  (0.45  cm/MW  for  InAs)  for  all  N 
until  approximately  10'^  cm"^  is  reached,  after  which  it  de¬ 
creases  like  a  Fermi  function  for  the  given  T  and  could  be 
fitted  very  well  to  a  functional  form  of  I3=a  +  be~‘^^^^’',  where 
n  is  carrier  density  in  the  unit  of  10“'^  cm"^.  Of  the  param¬ 
eters  {a,b,c),  we  found  that  only  c  depends  on  T.  Nq,  which 

g 

is  normally  assumed  to  be  constant,  is  found  to  change  by 
two  orders  of  magnitude  when  T  is  changed  from 
80  to  500  K. 

The  FCA  rate  evaluated  with  full  bands  is  found  to  be 
proportional  to  N,  in  agreement  with  that  assumed  in  the 

g 

literature.  However,  the  proportionality  constant  cr,  which  is 
the  FCA  cross  section,  is  found  to  be  strongly  T  dependent. 
The  contribution  to  the  FCA  comes  more  from  holes  than 
from  electrons.  The  transition  from  light-hole  to  heavy-hole 
band  does  not  require  phonons,  whereas  the  absorption  of 
light  by  electrons  is  possible  only  with  phonons.  The  T  de¬ 
pendence  of  phonon  and  electron  distribution  results  in 
strongly  T-dependent  FCA  cross  section. 

We  modified  the  previously  developed  full-band- 
structure  calculation  of  the  lifetimes  associated  with 
electron-hole  recombination  through  Auger  and  radiative 
processes  to  address  the  issue  of  electron  and  hole  quasi- 
Fermi  levels.  The  calculated  values  of  the  lifetimes  could  be 
fitted  accurately  to  a  function  form  t= 

All  T-dependent  parameters  {N,c ,G,p,cr)  are  fitted  ac¬ 
curately  to  a  fourth  order  polynomial  in  T,  for  easier  and 

f.  a  and  /3  are  the  one-  and  two-photon  absorption  coeffi-  more  efficient  use  in  obtaining  self-consistent  solutions  to 
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^  =  -aIp-pil-cT{No  +  N)Ip, 
dz 


dN 

dt 

iff 

dt 


dipr 

dz 


2hv 


N 


ct*{Nq  +  N)1p  ^(3 fp 


■cT{No  +  Nff 


(1) 


(2) 


(3) 


(4) 


The  equations  relate  the  pump  and  probe  intensities.  Ip  and 
/pj,  the  concentration  of  the  generated  charge  carriers  N,  and 
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rate  equations.  The  details  of  carrier  density,  coefficient,  and 
recombination  rate  variation  with  T  and  N  and  the  fitted  pa¬ 
rameters  will  be  given  elsewhere.  Here  we  present  only  the 
transmission  results  obtained  using  these  calculated  values 
and  parameters. 

The  calculated  values  of  f3{N,T),  cr{T),  and  Tj({N,T)  are 
substituted  in  Eqs.  (l)-(4)  and  solved  for  Ip^{z,r,t)  and 
N{z,r,t)  by  the  finite  difference  method  with  the  boundary 
conditions  M2,?-,0)  =  0, 

7’(z,r,0)  =  7’iatt,  and  /pj(0,r,r)=/[J,.e"^''^''pr)^,  where  rp  and 
are  the  pump  and  probe  radii.  The  beam  is  assumed  to  be 
radially  symmetric,  and  the  beam  energy  E  is  the  time  and 
r-integrated  value  of  Ip{Q,r,t).  Care  is  taken  to  yield  a  nu¬ 
merically  stable  solution  even  under  illumination  of  a  strong 
pump  laser  pulse.  For  a  given  incident  pump  irradiance,  we 
obtain  the  value  of  the  probe  transmission  as  a  function  of  r, 
t,  and  z  in  the  sample.  Then  the  transmission  intensity  is 
integrated  over  r  to  obtain  a  value  at  a  given  z-  The  trans¬ 
mission  intensity  evaluated  at  the  exit  surface  {z=L)  is  com¬ 
pared  with  experimental  transmission  curve.  Although  the 
calculated  curve  has  trends  similar  to  that  observed,  we  find 
that  inclusion  of  the  SRH  mechanism  is  needed  for  quantita¬ 
tive  comparison  with  experiment.  The  SRH  scattering  is  not 
intrinsic  to  the  material  and  depends  on  external  variables 
such  as  growth  temperature  and  pressure.  Our  detailed  cal¬ 
culations  carried  out  previously^^  indicate  a  nearly 
T-independent  SRH  lifetime  of  200  ns  in  2  X  lO'®  cm'^ 
doped  InAs.  Although  the  InAs  sample  studied  here  is  only 
10*^  cm'^  n  doped,  the  excess  carriers  generated  by  the 
pump  is  of  the  order  of  (1-2)  X  10'®  cm“^  and  justifies  the 
use  of  the  above  calculated  SRH  lifetimes.  We  further  as¬ 
sume  that  SRH  lifetime  decreases  linearly  with  excess  carrier 
density  and  obtain  /pr(L ,  t)  in  InAs  for  three  values  of  E.  The 
substrate  ambient  temperature  of  300  K  is  used. 

First,  the  equations  are  solved  only  for  transmitted  pump 
intensity  without  any  probe  present.  We  see  in  Fig.  1  that  the 
calculated  values  (solid  lines)  as  a  function  of  input  intensity 
agree  very  well  with  the  values  measured  for  four  different 
InAs  thicknesses.  Then,  for  a  fixed  InAs  thickness  and  three 
different  beam  energies,  the  transmitted  probe  intensity  is 
calculated  and  compared  with  experiment,  as  shown  in  Fig. 
2.  The  calculations  explain  the  observed  variations  very  well. 
It  is  important  to  note  that  good  agreement  is  obtained  for  the 


pump  transmission,  probe  transmission  recovery,  and  the 
minimum  of  the  probe  transmission,  which  depend  critically 
on  all  quantities:  the  intrinsic  carrier  density,  yS,  cr,  tempera¬ 
ture,  and  the  lifetimes.  Only  with  N  dependence  of  SRH, 
included  here  in  an  approximation,  an  excellent  agreement 
with  experiment  is  obtained.  The  slight  disagreement  at  large 
f  is  possibly  because  of  actual  N  dependence  of  SRH  life¬ 
times,  not  considered  here. 

In  conclusion,  we  have  used  pump-probe  experiments  to 
understand  the  optical  interaction  in  InAs.  The  results  are 
explained  with  accurately  calculated  TPA,  FCA,  AR,  and  RR 
lifetimes  and  iterative  solutions  to  coupled  rate  equations. 
Our  calculations  explain  the  observed  variation  in  output 
pump  intensity  and  time-  and  energy-dependent  probe  inten¬ 
sities  accurately.  We  showed  that  with  accurate  evaluation  of 
all  parameters,  the  modeling  can  provide  detailed  guidance 
in  understanding  the  underlying  mechanisms  and  in  choosing 
materials  for  light  propagation  applications. 
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